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Abstract 

We study localization-delocalization transition in quantum Hall systems with a random field of 
nuclear spins acting on two-dimensional (2d) electron spins via hyperfine contact (Fermi) interac- 
tion. We use Chalker-Coddington network model, which corresponds to the projection onto the 
lowest Landau level. The inhomogeneous nuclear polarization acts on the electrons as an additional 
confining potential, and, therefore, introduces additional parameter p (the probability to find a po- 
larized nucleus in the vicinity of a saddle point of random potential) responsible for the change 
from quantum to classical behavior. In this manner we obtain two critical exponents corresponding 
to quantum and classical percolation. We also study how the 2d extended state develops into the 
one-dimensional (Id) critical state. 

PACS numbers: 73.20.Fz, 73.20.Jc, 74.43.-f, 31.30.Gs, 76.60.Es 
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Celebrated quantum Hall effect (QHE) is realized in a 2d electron gas subiected to a 

HQ 

strong perpendicular magnetic field and a random potential [1|, 12] • The uniqueness of this 
phenomenon is in high precision of the plateaux in the Hall component and very rich physics 
in the interplateau transitions 31]. Here we will study the influence of the nuclear spin fields 
^ on the critical exponents in QHE. 

The rich physics of the random potential in quantum Hall systems could be roughly 
divided into spin independent and spin dependent (spintronics) electron scattering processes. 
Magnetic impurities perturb the QHE transport very strongly and will not be considered 
here. Recently sharply growing attention was attracted to the physics of the hyperfine 
interactions in the QHE. It was suggested theoretically ^| and observed experimentally 

m 

that the underlying nuclear spin structure can provide the microscopic information 
on the 2d electron wave functions and provide strong influence on the precision and other 
parameters of a QHE system. 

The interaction between electron and nuclear spins in heteroi unctions under QHE con- 

riQ 

ditions is due, usually, to the hyperfine Fermi contact interaction |7| . This interaction is 
represented by the Hamiltonian: 

Hint = -Inhli ■ He, (l) 

where 7„ is the nuclear gyromagnetic ratio, Jj is the nuclear spin and He is the magnetic 
field on the nuclear site, produced by electron orbital and spin magnetic moments: 

He = -gPY.^-^''^[re-R^. (2) 

e 

Here re is the electron radius- vector, Se is the electron spin operator, (3 = eh/rriQC is the 
Bohr magneton, g is the electronic (jf-factor and Ri is the nucleus radius-vector. 

It follows from Eqs. ([T]) and ([2]), that once the nuclear spins are polarized, i.e. if 
^S?;-^^ 7^ ' charge carriers spins feel the effective, time-dependent hyperfine field 
Bfif = B^j-exp {—t/Ti) (Ti is a nucleus rebcation time) which lifts the spin degeneracy even 
in the absence of external magnetic field. In GaAs/AlGaAs one may achieve the spin split- 
ting due to hyperfine field of the order of the Fermi energy 5|, iGj.The inhomogeneous nuclear 
polarization acts on the electrons as additional (to the scalar potential of the impurities) 
confining potential Vhj = —fXBBhj 

The nuclear spin polarization, once created, remains finite for macroscopically long times. 
Intensive experimental studies p, O] of this phenomenon in QHE systems have provided a 



more detailed knowledge on the hyperfine interaction between the nuclear and electron spins 
in heterojunctions and quantum wells. It was observed that the nuclear spin relaxation time 
is rather long (up to 10^ sec) and the hyperfine field acting on the charge carriers spins 
. ext.e,nely Ingh, up to lO^G &. The nuclea. .ela.at.on tune depends st.ongly on the 
vicinity to the impurity and its sign |3| . The presence of the impurity (long range potential) 
provides the necessary energy conservation in the spin-exciton creation process leading to 
the nuclear spin relaxation. We can, therefore, expect the following scenario: nuclear spins 
being polarized by some external field will then relax differently depending on whether they 
are close to maxima or minima of the scalar potential created by impurities. Therefore, they 
should affect strongly tunneling of electrons through saddle-point potential. 

When random potential varies smoothly (its correlation length is much larger than the 
magnetic length as, e.g., in GaAs heterostructures) a semiclasscial description becomes rel- 
evant: electrons move along the lines of constant potential. When two equipotential lines 
come close to each other (near a saddle point) tunneling is feasible. In this paper we inves- 
tigate how this picture will be affected by strong nuclear polarization. We find that scaling 
of the localiztion length is^modified (Eq. (8)), which is the main result of this work. 

In the network model ^9], electrons move along unidirectional links forming closed loops 
in analogy with semiclassical motion on contours of constant potential. Scattering between 
links is allowed at nodes in order to map tunneling through saddle point potentials. Prop- 
agation along links yields a random phase 0, thus links are presented by diagonal matrices 
with elements in the form exp(i0). Transfer matrix for one node relates a pair of incoming 
and outgoing amplitudes on the left to a corresponding pair on the right; it has the form 



In order for a system to be invariant, on average, under 90° rotation the transmission 
and refiection at the next neighbor node are interchanged, i.e. the transfer matrix has the 
same as in Eq. ([3]) form with a parameter e' = — e j^. In order to obtain this relation one 
simply interchanges and Z4 (see Fig. 1) and brings a new transfer matrix to the form of 
Eq. ( [3]) . We therefore describe scattering at the nodes indicated in Fig. 1 by circles with 
transfer matrix T(e) and at the nodes indicated by boxes with T(— e). 

The node parameter e is a relative distance between the electron energy and the barrier 
height. It is related to the physical quantities descibing the system 
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e={E-{n + ]^)E2-Vo)/E,, (4) 

where Ei measures the ratio between saddle-point paramters and magnetic field, E2 is a 
distance between Landau levels at strong magnetic fields, and Vq is a reference point of a 
scalar potential (see [l^ for details). It is easy to see that the most "quantum" case (equal 
probabilities to scatter to the left and to the right) is at e = 0, in fact numerical calculations 
9|] show that there is an extended state at that energy. Numerical simulations on the 
network model are performed in the following way: one studies system with fixed width M 
and periodic boundary conditions in the transverse direction. Multiplying transfer matrices 
for slices and then diagonalizing the resulting total transfer matrix , it is possible to extract 
the smallest Lyapunov exponent A (the eigenvalues of the transfer matrix are exp(AA^)). The 
localization length is proportional to 1/A. Repeating calculations for different system 
widths and different energies it is possible to show that the localization length satisfies 
a scaling relation 

In the QHE the thermodynamic localization length ^(e) ~ |e| and v = 2.5±0.5. This is the 
main result 9(] and it is in a good agreement with experimental data for spin-split resolved 
levels numerical simulations using other models 12| and semiclassical argument jisl. [l^ 
that predicts u = 7/3. 

It is possible to model classical percolation using CC model as well. It was shown pi 
that when the relative height of the barriers fiuctuate in the infinite range, the percolation 
becomes classical (no tunneling is allowed) and classical percloation exponent Ud = 4/3 is 
retrieved. On the other hand, when the fiuctuations are finite, their width acts as irrelevant 
p—flyQ and doe. not affect. 

In the present work we modify CC model in the following way. We expect that the 
presence of a polarized nucleus near a saddle point of the scalar potential will modify a 
tunneling parameter e in Eq. (3) by changing Vq to Vq ± Vhf. More, we also expect that 
due to different relaxation rates (in the vicinity of impurities of different signs) the following 
scenario can be realized: nuclei situated near different types of saddle point (nodes of the 
model) will be polarized in opposite directions, breaking, therefore, isotropy of the system. 
We model this situation by introducing a parameter < p < 1 describing the probability 
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FIG. 1: Network model with missing nodes 



that there is a polarized nucleus near particular saddle point. Due to the effect of high 
hyperfine fields descibed above, we, as a rough approximation, can expect that the barrier 
becomes "infinite", i.e. the transfer matrix at the node is now a unit matrix. On the model 
language it means that the quasiparticle stays on the same horizontal link (see Fig. 1), and 
isotropy of the model is therefore broken. Obviously, when p = 1 a 2d system is broken into 
M one-dimensional chains, and, due to the fact that there is no backscattering, all states are 
extended independent on energy e and system width M. We, therefore, expect the smallest 
Lyapunov exponent A = 0, in contradistinction to the "ordinary" 2d extended state, where 
A is finite, and infinite thermodynamic localization length is recovered only after finite size 
scaling. In this sense p = 1 case is close to a Id metal found for a dirty superconductors 
with broken time- reversal and spin- rotational symmetries jiol . 

Before we present numerical results, let us discuss the possible form for the scaling of 
the renormalized localization length. Now, when we have "wiped out" on average a fraction 
p of the nodes, a quasiparticle should travel larger distance (times 1/(1 — p)) in order to 
experience the same number of scattering events. Therefore, naively, one would expect that 
the effective system width is now M(l — p)~^ and the scaling is 



On the other hand, we should take into account that the "missing" node actually does not 




(6) 
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FIG. 2: Renormalized localizator! length at critical energy e = as function of the fraction of 
missing nodes p for different system widths. Solid line is the best fit 1.24(1 —p)~^'^. Dashed line 
is the fit with "naive" exponent v = 1. 

allow the quasiparticle to propagate in the transverse direction (we have chosen the system 
in such way that, if there is no scattering, the quasiparticle stays on the same horizntal 
link). Usually for CC model and its generalizations the typical value of the renormalized 
localization length for the extended state is of the order of 1, meaning that in the extended 
state the quasiparticle is able to traverse the system of the width M. Therefore, in the 
present situation we could expect even larger value of in the extended state, i.e. 
dependence with u > 1. 

In order to find both critical exponents we start by studying a p-dependence for e = 0, 
corresponding to the development of a 2d extended state into a Id extended state. The 
results for system widths M — 16, 32, 64 are presented on Fig. 2, allowing the following fit 

^ = (l-prV(0), (7) 

where /(O) is the value of the renormalized locahzation length in the extended state e = 
for the standard CC model {p — 0). This value for the critical exponent is suspiciously close 
to the classical percolation exponent — 4/3. We also show visibly worse fit of the data 
with the "naive" critical exponent i/ = 1. 

We next use the value z/ found in Fig. 2 and study numerically renormalized localization 
length for various e and p < 1. All our data collapse on one curve with abscissa in the 
form M/^(e) where thermodynamic localization length diverges as ^ ~ e~^« with quantum 
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FIG. 3: Data collapse for all energies e, system widths M and all fractions p 7^ 1 of missing nodes. 

percolation exponent Ug 2.5. The results of the scaling are presented on Fig. 3. 

We argue that one can understand the appearance of the classical percolation exponent 
in Eq. (7) by considering a quasiparticle on the standard CC model deeply into the localized 
regime. In this case localization length is M- independent, meaning that a quasiparticle 
does not "feel" the boundaries of the system, and its thermodynamic localization length 
C — ^M- Therefore, a quasiparticle travels on the perimeter of the classical cluster of the 
typical size ^. Then by increasing the fraction p of the missing nodes, we increase the size 
of the classical cluster, actually making it infinite as p approaches 1. Therefore, (1 — p) acts 
as energy in the classical percolation problem, explaining the value 1.3 ~ 4/3. 

Finally, all numerical data we have obtained supports the following scaling relation 



We stress that this is the first result for the network models to produce both quantum 
and classical percolation exponents form the same data. To summarize, we have studied 
the influence of nuclear spins on the localization-delocalization transition in quantum Hall 
systems. We have found that the fraction p of polarized nuclei acts as a relevant parameter, 
leading to a new scaling relation for the localization length (Eq. (8)). 

One of us (V. K.) appreciates valuable discussions with Alexander Mirlin,Yuval Gefen, 
Baruch Horovitz and Yshai Avishai. I.V. is gratefuU to Mike Pepper, Cavendish, for hospi- 
tality during the summer 2005, when part of this paper was accomplished. I.V. acknowledge 
the EuroMagNET of FP6, RII3-CT-2004-506-239. 

Note added. When this paper was already submitted my coauthor Professor Israel Vagner 
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has died after a lengthy and courageous battle with cancer. 
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